sce <- readRDS("data/sce_original_RPs_MPs.rds")
de_res_b <- readRDS("data/de_results_baseline_original.rds")
de_res_a <- readRDS("data/de_results_stimulated_original.rds")
features <- "state"
plotPbExprsMod(sce,
k = "all",
features = features,
assay = "exprs",
fun = "median",
color_by = "subgroup2",
facet_by = "antigen",
shape_by = "activation")
CATALYST::plotExprs(sce,
color_by = "subgroup2",
features = features,
assay = "exprs")
features <- "type"
plotPbExprsMod(sce,
k = "all",
features = features,
assay = "exprs",
fun = "median",
color_by = "subgroup2",
facet_by = "antigen",
shape_by = "activation")
CATALYST::plotExprs(sce,
color_by = "subgroup2",
features = features,
assay = "exprs")
plotDiffHeatmapCustom(x = sce,
y = de_res_b$t_test,
col_anno = c("subgroup", "patient_id"),
k = "all",
top_n = 20,
fdr = 0.05,
lfc = 1,
all = TRUE,
eff_r = de_res_b$effect_size,
sort_by = "padj",
assay = "exprs",
fun = "median",
normalize = TRUE
)
print_DE_top_table(sce, de_res_b, method = "t_test")
plotDiffHeatmapCustom(x = sce,
y = de_res_a$t_test,
col_anno = c("subgroup", "patient_id"),
k = "all",
top_n = 20,
fdr = 0.05,
lfc = 1,
all = TRUE,
eff_r = de_res_a$effect_size,
sort_by = "padj",
assay = "exprs",
fun = "median",
normalize = TRUE
)
print_DE_top_table(sce, de_res_a, method = "t_test")
plotDiffHeatmapCustom(x = sce,
y = de_res_b$CyEMD,
top_n = 20,
fdr = 0.05,
col_anno = c("subgroup", "patient_id"),
k = "all",
lfc = 1,
all = TRUE,
eff_r = de_res_b$effect_size,
sort_by = "padj",
assay = "exprs",
fun = "median",
normalize = TRUE
)
print_DE_top_table(sce, de_res_b, method = "CyEMD")
plotDiffHeatmapCustom(x = sce,
y = de_res_a$CyEMD,
col_anno = c("subgroup", "patient_id"),
k = "all",
top_n = 20,
fdr = 0.05,
lfc = 1,
all = TRUE,
eff_r = de_res_a$effect_size,
sort_by = "padj",
assay = "exprs",
fun = "median",
normalize = TRUE
)
print_DE_top_table(sce, de_res_a, method = "CyEMD")
createVennHeatmap(de_res_b)
createVennHeatmap(de_res_a)
sce <- readRDS("data/sce_CD42b_RPs_MPs.rds")
de_res_b <- readRDS("data/de_results_baseline_CD42b.rds")
de_res_a <- readRDS("data/de_results_stimulated_CD42b.rds")
features <- "state"
plotPbExprsMod(sce,
k = "all",
features = features,
assay = "exprs",
fun = "median",
color_by = "subgroup2",
facet_by = "antigen",
shape_by = "activation")
CATALYST::plotExprs(sce,
color_by = "subgroup2",
features = features,
assay = "exprs")
features <- "type"
plotPbExprsMod(sce,
k = "all",
features = features,
assay = "exprs",
fun = "median",
color_by = "subgroup2",
facet_by = "antigen",
shape_by = "activation")
CATALYST::plotExprs(sce,
color_by = "subgroup2",
features = features,
assay = "exprs")
plotDiffHeatmapCustom(x = sce,
y = de_res_b$t_test,
col_anno = c("subgroup", "patient_id"),
k = "all",
top_n = 20,
fdr = 0.05,
lfc = 1,
all = TRUE,
eff_r = de_res_b$effect_size,
sort_by = "padj",
assay = "exprs",
fun = "median",
normalize = TRUE
)
print_DE_top_table(sce, de_res_b, method = "t_test")
plotDiffHeatmapCustom(x = sce,
y = de_res_a$t_test,
col_anno = c("subgroup", "patient_id"),
k = "all",
top_n = 20,
fdr = 0.05,
lfc = 1,
all = TRUE,
eff_r = de_res_a$effect_size,
sort_by = "padj",
assay = "exprs",
fun = "median",
normalize = TRUE
)
print_DE_top_table(sce, de_res_a, method = "t_test")
plotDiffHeatmapCustom(x = sce,
y = de_res_b$CyEMD,
top_n = 20,
fdr = 0.05,
col_anno = c("subgroup", "patient_id"),
k = "all",
lfc = 1,
all = TRUE,
eff_r = de_res_b$effect_size,
sort_by = "padj",
assay = "exprs",
fun = "median",
normalize = TRUE
)
print_DE_top_table(sce, de_res_b, method = "CyEMD")
plotDiffHeatmapCustom(x = sce,
y = de_res_a$CyEMD,
col_anno = c("subgroup", "patient_id"),
k = "all",
top_n = 20,
fdr = 0.05,
lfc = 1,
all = TRUE,
eff_r = de_res_a$effect_size,
sort_by = "padj",
assay = "exprs",
fun = "median",
normalize = TRUE
)
print_DE_top_table(sce, de_res_a, method = "CyEMD")
createVennHeatmap(de_res_b)
createVennHeatmap(de_res_a)